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Abstract 

For a class of singular Sturm-Liouville equations on the unit interval 
with explicit singularity a(a + l)/x 2 , a £ N, we consider an inverse spectral 
problem. Our goal is the global parametrization of potentials by spectral 
data noted by A", and some norming constants noted by n a . For a — 
and 1, X a x K a was already known to be a global coordinate system on 
Lj|(0, 1). With the help of transformation operators, we extend this result 
to any non-negative integer a and give a description of isospectral sets. 

1 Introduction 

Inverse spectral problem for Schrodinger operator with radial potential H := 
—A + g(||Jf||) acting on the unit ball of M 3 leads by separation of variables 
(see [10], p. 160— 161) to consider a collection of singular differential operators 
Ho (q) j a G N acting on Ljjj(0, 1), denned by 

E a (q)y(x) := (-^ + + q(x)j y(x) = Xy(x), x G [0, 1], A e C, 

(1.1) 

with Dirichlet boundary conditions 

2/(0) = y(l) = 0. (1.2) 

Our purpose is the stability of the inverse spectral problem for H a . For this, 
we construct for each a € N, a standard map A a x n a for potentials q G £r(0, 1) 
with spectral data A a and some terminal velocities n a . 

Borg [1] and Levinson [7] proved that A x k° is one-to-one on Lj|(0, 1). 
Poschel and Trubowitz [9] improved this result obtaining A global 
real-analytic coordinate system on Lg(0, 1). Guillot and Ralston [6] extended 
their results to A 1 x k 1 . Next Carlson [3] (see also Zhornitskaya and Serov [12]) 
proved that for all real a > —1/2, A a x n a is one-to-one on Ljj|(0, 1). 

We complete these works, to any integer a, with the realization of A a x K a 
as a local real-analytic coordinate system on Ljjj(0, 1). Moreover, we obtain a 
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description for the set of potentials with same spectral data : these so-called 
isospectral sets are proved to be regular manifolds, expression for their tangent 
and normal spaces are given. The key point of the proof lies in the use of 
transformation operators: they help us to handle Bessel functions naturally 
underlying in this problem. 



2 The direct spectral problem. 

Part of the following properties are deduced or came from [6], [3]- [2] and [12]; 
nevertheless main structure is given by [9]. 

Define uj := A 1 / 2 (the square root determination is pointless because of the 
parity of functions using it) and for all n G N 

(2n+l)!! := 1 • 3 • • • (2n + 1). 

A fundamental system of solutions for (1.1) when q = is given by 

(2a +1)!! co a 

U i X >X) = 3a{^X), V{X,X) = ~ ( 2a+ 1)11 Vai^x), 

where j a and rj a are spherical Bessel functions (see section 5). This family 
is called fundamental since its wronskian W(u,v) is equal to 1. From their 
behavior near x = 0, u(x, A) is called the regular solution, it is analytic on 
[0, 1] x C; v(x, A) is called the singular solution, it is analytic on (0, 1] x C. 

Following Guillot and Ralston [6] but also Zhornitskaya and Serov [12], we 
construct solutions for (1.1) by a Picard's iteration method from u and v. 

Let if and ip be defined by 

<p(x,\,q) = ^2<Pk(x,X,q), ii(x,\,q) = ^^ fe (x, X,q) 

k>0 k>0 

with 

(fo(x,\,q) = u(x,X), 

r (2 1) 

fk+i(x, X,q) = I g(x,t,X)q(t)ip k (t,X,q)dt, k e N; 

(2.2) 

Q is called Green function and is given by 

Q{x, t, A) = v{x, X)u{t, A) - u(x, X)v(t, A), (x, t) e (0, 1) x (0, 1). 
Now follows the expected result: 
Lemma 2.1. Series defined in (2.1), respectively in (2.2) uniformly converge on 



ipo(x,X, V) = v(x,X), 
i> k+1 {x,X,q) =- / g{x,t,X)q(t)$ k {t,X,q)dt, keN. 



bounded sets of [0, 1] x C x L^(0, 1), respectively of (0, 1] x C x L 2 (0, 1) towards 



solutions of (1.1). Moreover, they satisfy the integral equations 

tp(x, A, q) = u(x, X) + f g(x,t,X)q(t)(p(t,X,q)dt, (2.3) 
Jo 

${x,X,q)=v(x,X)- [ g{x,t,X)q{t)ii{t,X,q)dt, (2.4) 
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and the estimates 



\p(x,X,q)\ < C e l Im "l*( r 
ip(x, 



a+l 



+ \ u \ x , 

,X,q)\ < c e li— Id— ) ri+i^V, 
with C uniform on bounded set of L^(0, 1). 



Proof. The proof is given for p, it is the same for ip. Estimates (5.4) for Bcssel 
functions gives 

x 



|u(x,A)| < Ce^ Im ^ x 
Iterative relation (2.1) leads to 



1 + \uj\x 



a+l 



(2.5) 



(pi(x,\,q) = / Q(x,t, X)q(t)u(t, X)dt, 
Jo 

which, combining (2.5) and the Green function estimate (5.6), is bounded by 

' t\q(t)\ 



\Mx,\q)\<C 2 e^ 



1 + \uj\x 



o 1 + M* 



dt. 



By successive iterations and recurrence, for all positive integer n, we get 



\(p n (x,X,q)\ < 



C 



n+1 



J Im ai\x 



1 + \uj\x 



a+l 



t\q(t)\ 

o i + M*' 



dt 



This proves uniform convergence on bounded sets of [0, 1] x C x £^(0, 1) for ip 



and the estimate. Integral relation follows from (2.1). 
Corollary 2.1. p> and ip, solutions for (1.1) follow the estimates 

; t\q(t)\ 



\p(x,X,q)-u(x,X)\<Ce^ x [jT\^\ 
$(x,X,q)-v(x,X) < C e^^ (1±ME 



o i + M* ' 
1 *k(*)l 



Le< ' denote the derivative with respect to x, we also obtain 

ip (x, X, q) = u (x, X) + J —{x,t,X)q(t)ip(t,X,q)dt, 

ip (x, A, q) = v (x, X) - / — {x,t,X)q(t)ip(t,X,q)dt, 



dt, 
dt. 



□ 

(2.6) 
(2.7) 



and 



\,'(x,X,q)-u'(x,X)\<Ce^^^) a £- 



t\q(t)\ 



dt 



X,q) - v'(x,X)\ < Ce)^-*) ^±^_J +1 £ _ 
with C uniform on bounded sets o/Lp(0, 1). 



Mi 
*k(*)l 

+ \uj\t 



(2.8) 
dt, (2.9) 
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Proof. For the first estimates, we just have to add those from the proof of the 
preceding lemma beginning at k = 1. Integral equations follow by derivation and 
the last estimates follow using estimates (5.7) and (5.9) for Q given in annexe 
5. □ 

Remark 1. Notice that for u^O, 

t 1 \uj\t 1 

< = — — < — 

~ 1 + \u\t M 1 + \oj\t ~ 

Thus, bound from the corollary are asymptotic estimates when u> — ► oo. 
Remark 2. Theses estimates lead to 

W(X,q)~W(<p(-,X,q)M-,X t q)) = l + o(^j, w -+ oo. (2.10) 

According to [8] and [2], following [6], there exists a function i/j(x,X,q) such 
that x a ip(x, X, q) is analytic on C x Lg(0, 1) for all x G [0, 1] and such that {ip, tp} 
is a basis for the solutions of (1.1), with the normalization 

W(9?,V) = 1- 

Behavior near x = for ip is thus inherited from u and with the wronskian 
relation this leads to the following boundary conditions at x = 



<fi(x,X,q) 

x->0+ X a+1 x^0+ 



lim y J a+ [ ' = 1, Um x a tp(x, A, q) = 1. 



Thus, ip is called as the singular solution of (3.1). 

With the estimate (2.10), W(A, q) has no zeros near the real axis and we can 
take for this singular solution the one defined by 

it \ ^ $( x iX,q) , 011 x 

^ X ^ = -W{x7W ( } 

For real values of a, it is interesting to read [3] and [12] (precisely when 
a > ^), particularly non integer values allow to study the spectral problem 
for the radial Schrodinger operator acting on the unit ball in a even dimension 
space (in our case, the operator acts on I 3 ). 

Uniform convergence for series in the previous lemma give regularity for tp 
and expressed as bellow 

Proposition 2.1 (Analyticity of solutions). 

(a) For all x £ [0, 1], y>(x, A, q) is analytic on C x L%,(0, 1). Moreover, it is real 
valued onlx £r(0, 1). 

(b) The map 

<p : (X,q) h-> <p(;\,q) 
is analytic from C x L^(0, 1) to # 2 ([0, 1],C). 

(c) For all x G (0, 1], ip(x, A, q) is analytic on C x L^(0, 1) and real valued on 

tx4(o,i). 
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Consequently, W(A, q) is an analytic function on C x L^(0, 1) and thus ip 
denned by (2.11) has the needed regularity. 

Regularity leads to existence of derivatives, precisely we obtain expression for 
Lc(0, l)-gradients with respect to A and q. For this, following [9], we calculate 
differential for <p with respect to q: 



[d q ip(x,X,q)](v) = / Q(x,t,X,q)v(t)(p(t,X,q)dt, 
Jo 

where 

Q(x, t, X, q) = ip(x, A, q)<p(t, A, q) - ip(x, X, q)ip(t, A, q). 
Rewriting this relation as a scalar product, we obtain 

Proposition 2.2. [Gradients for the regular solution] 
For all v € L^(0, 1), we have 

V q <p(x,X,q)(t) = ip(t,X,q) [ip(x,X,q)ip(t,X,q) - ip(x, A, q)ip(t, A, q)} l[ , x ](t), 

— (x, X,q) = — [d q (p(x, A, q)} (1), 

W q (f'(x,X,q)(t) = ip(t,X,q) [tp'(x,X,q)ip(t,X,q) - ip'(x, A, q)ip(t, X, q)] l[o,x](t)- 

2.1 Spectral data 

From this point, potential q has real values, that is 
2.1.1 Eigenvalues 

Spectrum localization and normalization for eigenfunctions is made following [ ] 
and [6]. Boundary conditions given by (1.2) define spectra of (1.1)-(1.2) as the 
zeros set of the entire function A i— > ip(l, A, q). Results from Zhornitskaya and 
Scrov [12] give simplicity and localization similarly to [9] and [6]. Regularity 
and expression for the gradient for each eigenvalue follows from [2] . 

Notations. Let {X a , n {q)) n>1 be the set of eigenvalues representing the spec- 
trum H a (q). For every integer n > 1, let g n (t, q) be the eigenvector with respect 
to X a , n (q) defined by the normalization 

9n(t,q) = 



«)ll£j(0,l) 

Now recall previous results with the following : 

Theorem 2.1. Let q G £r(0, 1), the spectrum for the problem (1.1) with Dirich- 
let boundary conditions (1.2) is a strictly increasing sequence of eigenvalues 
X a (q) = (A Qi „ (<?))„>! which are all real-analytic on L^(0, 1) and verify 



V 9 A 0)7l (i) = g n (t,q) 2 . 
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Corollary 2.2. 

Let q £ £r(0, 1), the following asymptotics 



g n (t, q) = V2j a (u a ,n(q) t) + 0\^A 



^ q^a,n 

(t) = 2j a (iu a , n (q)tf +0 

A«,n(?)= (n+^\ 2 + 0(l). 
are uniform on bounded sets in £j|(0, 1). 
Proof. Rewriting relation (2.6) gives for weR, 



(2a + 1)!! 
ip(x, A, q) = ^-^j^ux) + O j 



(2.13) 

(2.14) 
(2.15) 



u) \1 + \uj\x 



a+V 



Thus, 



1 , x , 2 , [(2a + 1)!!] 2 f 1 



,„ / 1 / \ui\x 
\ w \ 1 + |a>|x 



(2.16) 



dt. 



Estimate (5.4) for Bessel function j a leads to 

2 



a {wt) + Ol - ' 



=j„(wtr+0 - 



1 / wla: 



u) \1 + \uj\x 



2a+2 N 



then 



^(t,\,q) 2 dt 



[(2a +1)!!] 2 



j a (Ljtydt + o[ - 



Relation (5.14) gives 



1 , , , 2 , l[(2a+l)!!] 2 



1 + - 

UJ 



(2.17) 



Using both (2.16) and (2.17), knowing that u a ^ n {q) = \ n + — j tt + 0(1), rela- 
tions (2.13) and (2.14) follow. 
For (2.15), we write 



K,n(q) - A o , n (0) = 



1 d 



— (Xa,n(tq)) dt = j (VtqK,n,q) dt, 



[04] 2 



2ja(ua,n(tq)x) q(x)dx dt + 0[ - 



Thus, 



(q) = A a ,„(0) + / q(t)dt 



[0,1] 2 



2j a (u a , n (tq)xy ~ 1 q{x)dxdt + 0[ - 1 . (2.18) 
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Eq. (5.4) shows the boundedness for the above second integral. According to 
[11] Eq. (2.10), we have 

A.,«(0) = (n + |) - a{a + 1) + l 2 (n), 

which proves (2.15). □ 

The asymptotic for eigenvalues is far from optimal, but at this point we 
cannot do better unless more work to deal with the second integral in (2.18). 

2.1.2 Terminal velocities 

Solving inverse spectral problem implies the knowledge of other spectral data. 
Indeed even for the regular case, the eigenvalues don't give the associated poten- 
tial q. Nevertheless, partial results (as Borg theorem) exist. Thus, let's define 
the needed additional data as in [9] and [6]. 

Definition 2.1. Let K a , n , called terminal velocity, be defined for all integer 
n > 1 by 



K a ,n{q) = In 



u'(l,A a , n (0)) 
Moreover, let a n (t, q) be the function given by 

a„(t, q) = ip(t, x a ,n(q),q)i>(t, x a ,n(q), ?)• 

Give some properties for these items: 

Theorem 2.2. For all n £ N, q i— > n a ,n{q) is real- analytic on Lj|(0, 1). Its 
gradient is expressed by 

V 9 K a ,n(*) = -a n {t, q) + V q Xa.n(t) / a n (s,q)ds, (2.19) 

Jo 



and follows the estimate 



-ja{Ua,nt)?la{uJa,nt) + O ~r , 



(2.20) 



uniformly on bounded sets in [0, 1] x L§(0, 1). 



Preuve. Regularity of n a , n comes from regularity of A aj „ and ip. Expression for 
its gradient is a straightforward calculation using Proposition 2.2. 

Determine the asymptotic value of the gradient. Estimations (2.6)-(2.7) and 
(2.10) lead to 



a n (x,q) 



1 



1 + \U) a . n X 
v(x,X a ,n) + 



( 1 
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then, definitions of u and v give 



a n (x,q) 



\u a ,n\x 



1 + \(J a ,nX\ 
Tla{0Ja,nX) + O 



1 + \cj a , n x\ 



\u a< n\ x 



From (5.4) and (5.5), we deduce 



a n (t,q) = ja(u a ,nx)r] x)+0 



Now, we have 



a n (t, q)dt 



to, 



a, n J o 



t)Va(Wa, n t)dt + O 



Relation (5.15) gives the result. 



□ 



As noticed previously, on the contrary to the regular case, good estimations 
for terminal velocities are harder to obtain. Indeed, we have the relation 

Ko,n(<?) = J ^ («a,n (*<?)) & = J (V tq Ka,n,q) Lj(0,l) d *' 

Using Eq. (2.20), we get 

(tq)x)r) a (oj a ,n(tq)x) 



«a,n(?) = 

Thus, follows 



[0,l] 2 



U a ,n(tq) 



«a,n(<?) =0[- 



qix)dxdt + 0[ — 
n 2. 



(2.21) 



losing the C(^) accuracy. As Guillot and Ralston [6] for a = 1, we'll use 
transformation operator to avoid this problem. 

We finish our study of the direct problem some properties of the gradients. 

2.2 Orthogonality 

Proposition 2.3. For all (n, m) € N 2 , n,m>l, we have 



dx 



3- ( a„, -^; a m ) = 0. 
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Proof. Formally, the proof follows as in [9], thus we just prove the first one. Let 
(n, m) £ N 2 , with n, m > 1, integration by parts gives 



2 2 

ax 



1 



An 2 (ffm 2 )' - g m 2 (gn 2 )' dt= g n g m W(gn,g m )dt. 



When n = m. the first relation is true. Now suppose n ^ m, we use the identity 



T"VV(ff n ,5 m ) = (V« - K.m)g n gn 
ax 



to obtain 



2 2 
(IX 



□ 



Now, we deduce algebraic properties for gradients. For this, recall the defi- 
nition of the linear independence in an Hilbert space. 

Definition 2.2. An infinite vector family {uk)k>i m an Hilbert space is called 
free or its vectors are linearly independent if each vector of this family doesn't 
belong to the closed span of the others. Precisely 



Vfc > 1 , u k i Vect {uj \j >l,j^k). 
Corollary 2.3. [Gradients and orthogonality] 

(a) 1, {g n 2 — l} n >i are linearly independent, and so are .9" 2 |' 
over, these two families are orthogonal. 

(b) For all (n, m) S N 2 , n, m > 1 we have 

d_ 

dx 



More- 



(i) ^ VqK aj n, — — (VqK ai7n ) 



= o, 



(U) ( VqK aj „, — (V q \ a ,m) 



(Hi) ( V g A 0! „, — (VgA a>TO ) 



2*" 
0. 



2.3 The spectral map 

We define the spectral map A a x n a : 1^(0, 1 ) -> K x £j|P x with 

[A a x K a ] (g) = (J 9 (t)dfc, (Ao.nC?))^!, Kt?))^ 
where A 0l „(g l ) be given, through (2.18), by 



(2.22) 













I 









Regularity for A a x K a follows as in [9] (see also [G] when a = 1). We give this 
property with the following 
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Theorem 2.3. The map \ a x k" is real-analytic from L^(0, 1) inlx^X £j 
7fc differential is given by 



Keep in mind that this result is not optimal for the arrival space. To improve 
this, we use the following transformation operators. 

3 Transformation operators 

This idea is the key point for such an inverse spectral problem. It was first 
introduced by Guillot and Ralston for a = 1. Their goal was to transform scalar 
products with the first Bessel functions into scalar products with trigonometric 
functions. Successfully used in [2] and [4], Rundell and Sacks in [11] present 
theses operators stepwise for any integer a: First they construct elementary 
operator S a which "maps" Bessel function related to H a (q) into Bessel function 
related to H a _i(<j). Then, they just chain these operators to reach Bessel func- 
tion for a = 0, in other words, trigonometric functions. The following results 
extend their results. 

Notations. Let $ a and ^ a be defined by 



Our first add, as Guillot and Ralston for a = 1, is a transformation property 
for the family (^ a )a related to the terminal velocities. A second and new im- 
provement is some transformation properties for the families (& a ')a an d (^a')a 
related to the dual family appeared in corollary 2.3. We use it further to give 
precision for isospectral sets. 

Lemma 3.1. For all a £ N, a > 1 define S a acting from L^(Q, 1) in L^(0, 1) 




®a{x) = ja{x) 2 et ® a (x) = j a (x)r] a (x), ire [0,1 



by 




(3.1) 



We have the following properties: 
(i) The adjoint operator for S a is 




(3.2) 



(ii) The family {S a } pairwise commute: 



V(a,6)eN 2 , S a S b = S b S a . 



(Hi) S a is a bounded operator in 



■l(o,i). 




Its inverse is given by the bounded operator on 



(0, 1) defined by 
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(v) $ a and ^ a check the properties 

$ a = -S* [$ B _i], tf a = [*„_!], (3.3) 

*i=-^[tl], < = -^a[*a-J- (3-4) 

Lemma 3.2. For all a e N* Zei T a 6e defined by 

T a = (—l) a+ S a S a -i ■ ■ ■ S±. 

Then 

(i) T a is a bounded one-to-one operator on L^(0, 1) such that : 
For all q £ L^(0, 1) and all AeC, 



-- sin(2A:r) 



(3.5) 



(3.7) 



/ [2$ a (\t) - l]q(t)dt = f cos(2\t)T a [q]{t)dt, 
lo Jo 

1 f 1 

* a {Xt)q(t)dt = -- / sm(2\t)T a [q](t)dt. (3.6) 
* Jo 

(ii) The adjoint of T a check 

2$ a (Az) - 1 = T* a [cos(2Az)] , * a (Ax) = T a * 
and its kernel is 

ker(T a *) = Span {x 2 , x 4 , . . . , x 2a ) . 

(Hi) T a define a Banach isomorphism between L^(0, 1) and (ker(T*)) . 
/is inverse is given by the bounded operator in L^(0, 1) defined by 

B a [f] := {-lf+ 1 A a A a _ l ---A l , 

moreover, we have 

<(A.t) = B a [- sm(2Ax)] , & a (\x) = B a [- cos(2A.t)] (3.8) 

Proof of theses results lies on properties of Bessel functions for the calculation 
part and lies on Hardy inequalities for the bounded properties and range of 
presented operators. For a detailed proof, we send back to [11], the new facts 
follow similarly. 

The first use of these operators is to improve spectral data estimations. 
Proposition 3.1. Uniformly on bounded sets in Lj|(0, 1), we have 









(n-\ 




L 









nKa,n(q) = £ 2 (n). 
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Proof. According to (2.18) and (2.21), we just have to find the behavior for 
integral terms /„ and J n defined by 



In(q) 



[0,1] : 



2ja(ua,n(tq)x) -lq(x)dxdt, 



[0,1] 



J n (q) ■- 

We use (3.5) and (3.6) to get 
In(q) 



(tq)x)i] a Uo a ^{tq)x) 



q{x)dx dt. 



[0,1] 2 



u a ,n(tq) 

cos (2L0 a , n (tq)x)T a [q](x)dx dt, 



Uq) = - I! 

J J\0,l] 



sin (2uj a , n (tq)x) 



'[o,i] 2 ^a,n{tq) 
From (2.12), we have u a , n = (n + § ) tt + O(^). Hence 



T a [q] (x)dxdt. 



I„(q)= J cos ((2n + a)irx)T a [q](x)dx +o(J-^ 



Jn{q) 



1 



in ((2n + a)nx) T a [q](x)dx + O — 



(2n + a)n J 
These relations lead to the result. 

Give two more estimates for the following functions 

Notations. For all integer n > 1, we define 



□ 



V ajn {x,q) := 2 



d_ 

dx 



^ q^a,n 



W a ,„{x,q) := -2 



d 

dx 



V qf^a : n 



(3.9) 



Proposition 3.2. Uniformly the bounded sets in [0, 1] x L^(Q, 1), we have 



V a ,n(x, q) = ^UJa.n (®a'(u a , n x) + O f J 

W a , n (x, q) = -2^ a '(uj a . n x) + o(^j . 



(3.10) 
(3.11) 



Proof. We only proof the first one, the second follows similarly. Wc have, ac- 
cording to the eigenvalue gradient expression, 

V a ,n(x,q) = 2-^j- (g n (x,q) 2 ) = 4g n (x,q)g' n (x,q). 
Estimate (2.8) combined with (2.17) and (2.10) leads to 

g n '{x, q) = V2~U>a,nja{Ua,nX) + 0(1) . 

We conclude that to (2.13) which allows us to write 



V atn (x,q) = 4iLl a . n ( 2j a (LO a , n x)j a ' (<jj a , n x) +0[- 



□ 
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4 The inverse spectral problem 

Now, the map A° x K a defined by (2.22) is a real analytic map from L^(0, 1) 
into R x £g x ^jj. Give the main result of this paper: 

Theorem 4.1. For a/Z integer a > 1, d q (X a x K a ) is a Banach isomorphism 
from Lj|(0, 1) onto R x £g x 

From the results of Zhornitskaya, Zerov [12] and Carlson [3], for all a > —i, 
the map A" x k° is one-to-one on Lj|(0, 1). We immediately deduce the following 
property 

Corollary 4.1. For all integer a > 1, i/ie map A a x n a is a global real analytic 
coordinate system on L™(0, 1). 

Proof of the theorem. The main arguments of this proof are extensions of those 
in [6]. From orthogonality relations in corollary 2.3 the family 

{i}u(v g A ,„) U {V 9 K a ,„} n>1 

l J n>l — 

is free in £j|(0, 1). Let r n and s n be defined by 

r n (x) = V q \ a ,n( X ) ~ (2$a(u a , n x) ~ 1) , (4-1) 

s n (x) = V q n a . n (x) — *k a {uj a>n x). (4.2) 

Using lemma 3.2, we have, for all v £ Lj|(0, 1), 

(V q \ a , n) v}= f {cos(2Lu a , n t) + R n (t))T a [v}(t)dt, (4.3) 



-1 



(V g Ka,n,v) = - / (sin (2w , n t) + S n (t)) T a [v](t)dt, (4.4) 

^LOa,n Jo 

with 

R n = B* a [r n ]et^=B* a [s n }. (4.5) 
Moreover, recall that T a [l] = — 1, particularly 

(l,v) = [ v{t)dt= [ -T a [v](t)dt. 
Jo Jo 

We denote F the operator defined by 



F(w) = ^ {-l,w), {(cos (2cj a , n t)+R n {t),w)} n ^ 1 



—n 



(sin (2w , n t) + S^C*)) , w 



such that d 9 (A a x n a )(v) = F o T a [w]. From lemma 3.2, T a is invertible from 
.Lj|(0, 1) onto (Span {a; 2 , x A 1 . . . ,x 2a }) ± . Thus, we have to show that F is in- 



vertible from (Span {x 2 , x 4 , . . . , x 2a }) onto R x ^ x Following [6], we 
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have to show the invertibility of F the operator from L^(Q, 1) into Ixl^x 
sending functions in L^(0, 1) on their Fourier coefficients (or, more simply, their 
scalar products against each element of the considered family) with respect to 
the family 

T = {1} U {i 2i } iG [i, o] u i cos ( 2w «.™*) + (*)>»>! 



U <j (sin (2w a , n t) + S n (t)) |> . (4.6) 

,n 



n>l 



This last statement will come from [9]: Appendix D, theorem 3: 

Lemma 4.1. Let T = {/„} ne n be a sequence in an Hilbert space H , with the 
two following properties 

(a) there exists an orthogonal basis £ = {e„} nS N in H for which 



Eli/-- 



e„|| 2 < oo. 



(b) the f n are linearly independent. 

Then {/n} rag N * s a basis of -ff and the map F : a: h- > {(/„, x)}„ 6 n is a linear 
isomorphism between H and £j|. 

Considering the orthonormal basis 

COS 7TX, sin7ra;, . . . , \/2cos (2n + l)7rx, \/2sin (2n + l)7ra, . . .| , 

when a = 2a + 1 with oeN (as in [6] for a = 1), and 

£ = |l, \/2 cos 27T2;, \/2sin27rx, . . . , y/2 cos 2n7rx-, \/2 sin 2n7ra;, . . . j , 

when a = 2a with a £ N (as in [9] for a = 0), estimates (2.14) and (2.20) used 
with (4.1) and (4.2) show that 

then the boundedncss of B* and relations given in (4.5) give 

ll^!lL i( 0,l)= O (^)> 11^(0,1) =°(£ 

which implies condition (a) after normalization. Lemma 4.2 will give condition 
(b), thus the proof will be complete. □ 

Now, give the lemma implying the theorem. Ingredients of the proof are 
similar to [6], we improve the method in a systematic way (for each a). 

Lemma 4.2. T defined by (4.6) is a free family in L^(0, 1). 
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Proof. Rewrite relations (4.3) and (4.4) as 

T* (cos (2w«, n t) + R n {t)) = V 9 A a ,„ 

and 

T a \7T~~ ( sin (2w tt , n t) + SnW) ] = V g K a , n . 

Since T* is bounded and the famille {1} U < V ' q \ a ,n \ U {V g K ai „} n>1 is free, 

I ' J n>l ' — 

we deduce the linear independency for the following family 

{1} U {cos (2w , n t) + i?n} n>1 U { — — (sin (2w„, n t) + S n ) 



2u> a ,n 



n>l 



Let k G [l,a], denote Wk the function Wfc(t) = t 2k . Let us show that Wk is 
not in the closure of Vect [T \ {Wfc}). (We might show iteratively that Wk £ 
Vect(J r \ {Wj,j G [fc, a]}), but it is not necessary by taking dm = for any 
to G [fc, a] in the next expression.) 

Suppose the contrary: there exists for j € N a sequence of vector 

W k U) (t) = a { j) + °$W m {t)+ Y, alp (cos(2uj a , n t) + R n (t)) 

m£[l,a] nG[l,JVj] 

+ E (sin (2a, a , n t) +£„(*)), 

neH,JVi] ^ a <" 

with TV,- < oo, a£\ b& G K such that 

wf } ^ in 4(0,1). 

Since T*(W m ) = for to = 1, . . . , a, the sequence 

w« := T* a {Wi j) ) = -a W + Y ^ViKn + bPV 9 Ka, n 

tends to in Lj|(0, 1) when j — > oo. Thus, corollary 2.3 gives 

atf' = / w^\t)dt — ► 0, (4.7) 
Jo 

alf = -2 [ w& (*)-£■ (V q K a , n ) dt — > 0, (4.8) 

- -2 f (V 9 A a ,„) A — ► 0. (4.9) 

Jo J-* 00 

Now consider w £ Cq°([0, 1],R), supported in [5, 1] with 5 > 0, such that 
(w, W m ) = 4,m, m G [1, a] 

and 

(B a [w],l) = i.e. (w,l) = c(A;,a)(a;,W fe ). 1 



'c^a) represents the Wk component of It can by computed by induction on a: 

B*[l] = -1 + Em=l c(m,a)W m - Particularly c(l,a) = a(a + 1) (idem [6] when a = 1). 
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Smoothness and support of ui imply that 



Lo(t) cos {2u a , n t)dt, / uj{t) sin {2uj a . n t)dt = — , ViV G N 



and that B a [u} is C°°([0, 1],R) supported in [S, 1]. 

Now, plug estimation (2.6) in the integral expression (2.3), then use (2.17), 
controls (5.4) and (5.6) to obtain the uniform estimate on [0, 1] 



r n (x) 



2u(x,X atn ) I Q(x,t,\ a , n )q(t)u(t,\ a , n )dt 
a+1 \ 



1 



2(Aa : „) 



u(x,X a ,n) +° 



(4.10) 



Vim;! (( 2 «+i)!0 2 

Thanks to lemma 4.3, this implies that the following family is summable 
{(w, [cos {2cj a . n t) + R n (t)])} n ^. 
Now turn to s n . From relations (4.2) and (2.19), we get 



s n {x) = -a n (x,q) + V q X a ,n(x) / a n (t,q)dt ^ a (uj an x), 

Ua.n 

-a n (x, q) + (2$ a (w aj „x) + r n (x)) j a„(t, q)dt — ^ a {u a , n x), 



= -a n (x,q) + 2$ a (cj ai „x) / a n (t,q)dt 

Jo 

f 1 1 

+ r n (x) / a n (t,q)dt ^ a (u) ain x). 

JO ^a,n 

Respectfully insert (2.6) and (2.7) in integral expressions (2.3) and (2.4); then 
use (5.4), (5.5), (5.6), (5.8) and wronskian estimate (2.10) to obtain the uniform 
estimate on [0,1] 



S n {x) = - v(x, X a ,n) / G(x,t, A tt , n )g(*)u(i, X a ,n)dt 



+ u(x,X a ,n) / G(x,t, X a ,n)q(t)v{t, \ a ,n)dt 



x) I a n (t,q)dt 



> S n (x) 



+ r n (x) / a n (t,q)dt 
Jo 

+ — * a ( Wa , nI ) (W^ 1 - 1) + O (\ ) . 
With the help of lemma 4.4, we deduce the summability of the family 

r 1 -i 

n oj(t)- (sin(2w Q ,„i) + S n (t))dt 

Jo 2w a n 



(4.11) 
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We finish the proof writing 

(w,Wjp}=ag ) (cj,l)+ ]T oW ) <w,(cos(2w , n t) + iJ n (t))) 

«e[i,w 3 ] 



2u a 

ne[l,JV>] 

Indeed, the following estimates, deduced from (3.10)-(3.11) and (4.7)-(4.9), 
\a$\ < C, |6W)| < Cn, V(n,j) G N 2 , 

imply that 

Uw®) o, 

which is contradictory with the choice of uj. □ 
Lemma 4.3. 

Proof. Let r„ i and r n: 2 be the first and second terms in(4.10). Recall that 

(uj,R n ) = (Lo,B* a [r n }) = {B a [w],r n ), 

thus, we just have to show that {(B a [uj],r n j)} € j = 1,2. First consider 
r„.2: 



(B a [o;],r n , 2 ) = ( B a [u], 



1 2(A„,„) 



a+l 



^a+l)!!) 2 )^'^ /' 



((2a +1)!!) 2 



B aH, — ,„+!,'' „ 2 - 1 2j a (w aj „x) 



2 



2(A a ,„) a+1 ||^ 



n II 2 



With (2.17), we have 

(B a [u],r n>2 ) =o(—) (B a [u],2 3a {Lo a . n xf) , 
then, hypothesis upon -B a [w] gives 

(B a [uj},r n , 2 ) =o(—) (B a [u)],2j a (uj an x) 2 - l) . 



Lemma 3.2 leads to the result. 

Now consider r n i and write as before 



r n ,i{x) = 2 ( 2u(x, \ a , n ) I G(x,t,\ a ,n)q(t)u(t,\ a , n )dt 

\\<Pn\\o V JO J 

^3a(u>a,nX) / Q{x,t,\ ajn )q{t)ja,(w a>n t)dt. 



Wn\\ 2 

2((2a+l)!!) 2 
K„) 2o + 2 ||<^ 
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With (5.4) and (5.6) we get 



Q(x,t, X atn )q(t)ja(^a,nt)dt = O 
\^a,n 



and, thanks to (2.17), we have 



From the expression of G(x, i, A), we can write 
4 r 1 r x 

/ B a [uj](x)j X ) I q(t)ja(Wa,nt)r]a(uJa,nt)dtdx 

U a ,n Jo Jo 

/ B a [uj}(x)j a (uj a ^x)r] a (u} a ^x) j q(t)j a (<J a ,nt) dtdx 

w «,n Jo Jo 
1 



o 



11- 



But, £?aM is supported in [5, 1], 6 > 0. Relations (5.2) and (5.3) thus give 
(B a [v],r ni x) = 



B a [uj](x)(l - cos(2w ai „cc - an)) / q(t)j a (uj a . n t)-q a {u) a ^ n i)dtdx 

J a,n JO Jo 

f X 

B a [oj](x) sin(2w aj „a; - air) I q(t)j a {w a ^ n t) 2 dtdx 

W a ,n Jo Jo 

Integrating by parts in all terms having factors of cos(w aj „a;) or sin(w a .„x), we 
get 



(B a [u], r n i) 



B a [uj](x) / q(t)ja{u> a ,nt)Tla(Va,nt)dtdx 
J a,n JO Jo 



o\- 2 

n z 



Finally, interchanging the order of integration and using properties of the trans- 
formation operator, we obtain the concluding relation 



(B a [u>], r„ ; i) 



Lemma 4.4. 



-1 



'Ja,n Jo 



sin(2uj a . n t)T a 



1 1— > q(t) j B a [u>](x)dx 



□ 



{{w,S n )} n>1 ell 
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Proof. Recall 



(UJ,S„) = 2uJ a , n {uJ 1 B* a [s n \) = 2(jJa,n{B a [u],S n ). 

The proof is similar to lemma 4.3 for s n (see (4.11)): we change in the scalar 
product (B a [u)],s n ), 2$ a (A Q; „a;) by 2$ a (A ai „x) — 1 and use the transformation 
operator; for the two following terms in s n , as previously, we have to use the 
support of B a [u>], asymptotics deduced from (5.2) and (5.3), integrate by parts 
terms with cos (2oj a ^ n x) or sin (2uj a ^ n x) and finally invert integration order to 
use transformation operators. 

The term following s n in (4.11) is controlled thanks to r n \ denote s n the 
remaining term. But, we have 



{B a [w],S„) 



O -5 / S„H(t)j (A (K,nt) dt 



= o 



1 



B a [u)](t) sin {2u> a ,nt)dt + O 



1 



Thus < n (-B a [a>], s n ) > is in ^ and the proof is completed. 

I J n>l 



□ 



Now comes the corollary (see (3.9) for the definition of V n and W n ). 

Corollary 4.2. A a x K a is a local real analytic isomorphism on L^(0, 1). More- 
over, the inverse of d q {\ a x n a ) is the bounded linear map from lx/gx!| onto 
£j|(0, 1) given by 



n>l 



n>l 



{d q {\ a x K, a ))- 1 {i 10 ,ri,^ = m + y>nW a ,„ + V —v a 



Proof. The first part of the corollary comes from the theorem 4.1 and from the 
inverse function theorem. Now consider (?7o, ?7, £) G K x x £\ and define 



Since B a is bounded, estimations (3.10) and (3.11), with relations (3.8) give 



and 



-V a . n {x,q) = B a 



W atn (x,q) = B a 



sin(2aj a: „a;) + O 



4w Q ■ 



-2 COs{2iJa,nX) +0\- 







< 






n 


(n 


)] 







Definition of £ and 77, estimation of eigenvalues and boundedness of B a imply 
the convergence in Ljj|(0, 1) for the series defining u. 
From corollary 2.3, we have 

(l,ti) = rjo, 

and for all integer n > 1 

^V g A a ,„,u^ = 7i„, {nS/ q n a .n,u) = £„. 
Thus, we have c? g (A a x n a ) (it) = (7/0, ? ?, 0: which proofs the corollary. □ 
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We finish with the description of isospectral sets. For go G ^r(0, 1), we define 
the set of potentials with same Dirichlet spectrum as go, called isospectral set, 
by Iso(<7o,a) = {<? 6 ^r(0, 1) : A a (g) = A a (g )} ■ The new fact of the following 
result is to explicit tangent and normal spaces. 

Theorem 4.2. Let q Q G Ljjj(0, 1), t/ien 

(a) lso(qo,a) is a real- analytic manifold o/L^ (0,1) of infinite dimension and 
codimension, lying in the hyperplane of all functions with mean qo(t)dt. 

(b) At every point q in lso(qo,a), the tangent space is 

T q Iso(q ,a)= lj2— V *,n-£d 

\n>l U 

and the normal space is 

N q Iso(g , a) = li]o + Y^ Vn{g n 2 - 1) : (m, v) G R x 4 



Proof. It is straightforward from [9] . 

(a) The first part of the assertion comes from [2] Theorem 1.3, the second is 
direct from (3.1). 

(b) Since 

T 9 Iso(g ,a) - (d q (\ a X O)" 1 ({0 Rx ^} x 4), 

corollary 4.2 gives expression of the tangent space. According to corollary 
2.3, {l,g n 2 — 1 : n > l} is free, orthogonal to the free family (V^, n ) nG z- 
Thus, we have 



Vo 



Vn(g 2 n - 1) : (Vo, V) G R X 4 I C N q Iso(g ). 



TI>1 



Moreover, any vector orthogonal to {l,g n 2 — 1 : n > l} is, with regards to 
the Frechet derivative of A a , in the null space of d q X a . Thus, we get the 
other inclusion and then the proof. 

□ 

To finish, we recall that the characterization of the spectra of each operator 
H a (g) was obtained by Carlson ([2] Theorem 1.1). 

5 Bessel functions 

Spherical Bessel functions j a and r\ a arc defined through 

/ IT Z I 7T Z 

3a(z) = J yJa+l/aC*). Va(z) = (~1)V y J-a-1/2^), (5-1) 
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where J v is the first kind Bessel Bessel function of order v (see [5] for precisions). 
From [5] formulas (1) and (2) section 7.11 p. 78, they behave like 



3a (z)=sm(z-—)+0\-r-r-\, \z\ 



CO, 



*,(*) = cob (*-f)+o( {z 

The following estimates can be found in [3] 
• Uniform estimates on C: 



Im z\ 



\z\ — > 00. 



(5.2) 
(5.3) 



1 + M 

M 

• Estimations for the Green function G(x, t, A) when < t < x: 

a+l / 1 i i , \ a 
X \ / 1 + \Ul\t \ 



(5.4) 
(5.5) 



|S(x,t,A)| <C 



1 + |u>|a; 



exp (| Imcj|(a: — t)) , (5.6) 



^(x,t,A) 



< C 



a /1 + Mi\ a 



exp(|Imw|(x-t))- (5.7) 



l + \u\xj V * 
• Estimations for the Green function G(x, t, A) when < x < t < 1: 



|£(x,t,A)| <C 



1 + Mx\ a / i 



.c 



1 + \u\t 



a+l 



exp(|Imw|(i-a:)), (5.8) 



< c 



1 + \lu\x 



a+l 



i + H* 



a+l 



exp (| Imcj|(t — x)) 



Rewriting relations (54) — (56) in [5] section 7.2.8 pp.11 — 12, we get 

Xja'(x) = Xj a -l(x) - aj a (x), 
Xja-l'(x) = aj a -l(x) - Xj a (x), 

xi-]' a {x) = xi-] a -i{x) - ar) a (x), 
xr^_i(x) = ar) a -i(x) - xrj a (x). 



j a (vt) dt = 



1 







1 + 



ja(uft)lh{ut)dt = 0[ - 

o V w 



(5.9) 

(5.10) 
(5.11) 
(5.12) 
(5.13) 



We also deduce the following uniform estimates with w € M, \u>\ — > +oo 

r-l 



(5.14) 
(5.15) 
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